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Abstract 

We consider an electro-weak scale model for Dark Matter (DM) and radiative 
neutrino mass generation. Despite the leptophilic nature of DM with no direct 
couplings to quarks and gluons, scattering with nuclei is induced at the 1-loop 
level through photon exchange. Effectively, there are charge-charge, dipole-charge 
and dipole-dipole interactions. We investigate the parameter space consistent with 
constraints from neutrino masses and mixing, charged lepton-flavour violation, per- 
turbativity, and the thermal production of the correct DM abundance, and calculate 
the expected event rate in DM direct detection experiments. We show that current 
data from XENON100 start to constrain certain regions of the allowed parameter 
space, whereas future data from XENON1T has the potential to significantly probe 
the model. 
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1 Introduction 



The Standard Model (SM) is very successful in describing the fundamental particles 
of our world. The only solid evidence for its failure so far is the fact that neutrinos 
have mass, which is a necessity due to the observation of neutrino oscillations [1-4]. 
Furthermore, the standard cosmological model, the ACDM model, provides an excellent 
description of our Universe, with the exception that within the SM there is no viable 
candidate for a Dark Matter (DM) particle, which is an important ingredient of the 
ACDM model, supported by observations such as the rotation curves of spiral galaxies [5], 
WMAP CMB measurements [6] and gravitational lensing [7]. Hence, both neutrinos, as 
well as DM require an extension of the SM. Often these two phenomena are considered 
separately, since they might be manifestations of physics from vastly different energy 
scales. Here we adopt the hypothesis that neutrino mass and DM are related, and both 
emerge from physics at the TeV scale. In this respect models which generate neutrino 
masses radiatively [8-14] are intriguing. 1 Loop suppression factors and several powers of 
Yukawa couplings can bring the scale of neutrino mass generation down to the TeV, and 
symmetries required to stabilize DM may play a role for neutrinos, for example forbid 
tree-level mass terms. Recent works in this context can be found in refs. [18-39]. 

The so-called WIMP hypothesis suggests that DM interacts sufficiently with SM par- 
ticles in order to generate the relevant abundance due to thermal freeze-out from the 
primordial plasma. This motivates the direct search for DM in our galaxy by look- 
ing for the scattering of DM particles with nuclei in underground detectors. Several 
direct detection experiments are pursuing such searches, for example the CDMS II [40], 
XENON100 [41,42], CoGeNT [43], DAM A/LIBRA [44], CRESST-II [45], ZEPLIN-III [46] 
and KIMS [47, 48] experiments. In typical WIMP models DM interacts directly with 
quarks, providing DM-nucleus scattering at tree-level [49,50]. Here we are interested in 
so-called "leptophilic" models, where DM couples directly only to leptons, see e.g. [51]. 
Even in that case, DM-nucleus interactions can be induced at loop-level due to the ex- 
change of the photon [52]. The resulting effective interactions have been investigated in 
refs. [52,53]. In the following we will consider a model where the corresponding loop- 
diagrams induce a magnetic and/or electric dipole moment interaction [54-61]. 

We consider a model proposed by Ma [10], in which neutrino masses are generated 
through 1-loop interactions and the particles which propagate in the loop can be DM 
candidates, being leptophilic by construction. The DM phenomenology of the model 
and extended versions thereof has been studied in refs. [19-23,33,34] and prospects for 
collider searches have been studied in refs. [23,27]. We consider the situation that the 
lightest right handed neutrino is the DM candidate and the second lightest right handed 
neutrino is almost degenerated with the DM candidate. Under this situation, elastic DM- 
nucleus scattering is extremely suppressed and inelastic scattering induced by a lepton- 

lr This is only a small collection of references. There are many more models which generate neutrino 
masses radiatively, also including supersymmetry, see e.g. [15-17]. 
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loop coupled to the photon gives the dominant contribution to the event rate in direct 
detection experiments. We calculate the event rate in the model and compare it with 
XENON100, KIMS and DAMA data. The paper is organized as follows. In Section 2, 
we shortly review the model from ref. [10]. We discuss the constraints from neutrino 
oscillation data, lepton-flavour violation and the thermal production of the DM relic 
abundance. In Section 3, we discuss the inelastic scattering cross section in an effective 
theory approach and calculate the event rate. Moreover, monochromatic photons from 
the decay of the excited DM state are also discussed. We summarize and conclude in 
Section 4. Explicit functions needed in the effective theory approach are listed in the 
Appendix A. 



2.1 Neutrino masses and mixing 

The model proposed by Ma in ref. [10] is a simple extension of the SM, which correlates 
neutrino physics and the existence of DM. The added particles to the SM are three right 
handed neutrinos Ni (i = 1, 2, 3) and one inert Higgs doublet 77. In addition, a discrete 
Z2 symmetry is imposed: odd for the new particles and even for SM particles. The new 
invariant Lagrangian is 



where <fi is the SM Higgs doublet. The vacuum expectation value (VEV) of 77 is assumed 
to be zero, so that the discrete Z2 symmetry which guarantees the stability of DM is 
an exact symmetry. Thus Dirac neutrino masses are not generated through the Yukawa 
couplings in Eq. (1). After electroweak symmetry breaking, the SM Higgs <fi obtains the 
VEV (4>°) and Majorana neutrino masses are generated radiatively with the effective mass 
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C N = N^PrN, + (D^ {D^rj) - -JiNfPnNi + h ai £ a rf P R N + h.c, (1) 



and the scalar potential V(<p, 77) is 



VO, 77) = mj^ty + mffirj + y ($<j)) 2 + y (77*77) 2 

+A 3 (0 f 0) Wv) + A 4 {fa) faV) + y (<P j v) 2 + h.c, 



(2) 




(3) 



where Mi are the masses of the right-handed neutrinos Ni, M 2 
and the loop function I(x) is defined as 



ml + (A 3 + A 4 ) (<f> ) 2 




(4) 
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These relations hold for small coupling A5, which is needed in order to obtain the correct 
neutrino masses, see below. This assumption is justified since an extra U(l) symmetry 
appears in the limit of A5 — >■ 0. 

As shown in ref. [21], the close to tri-bimaximal mixing of the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) mixing matrix is achieved by adopting the following flavour 
structure for the Yukawa couplings h ai (rows are labeled by a = e, /i, r and columns by 
i = 1,2,3): 



hai = \ h h 2 h 3 \ . (5) 




This matrix implies #23 — 7r /4, #13 = and tan# 12 = -^h' 3 /h 3 . 2 From the current best 
fit value sin 2 fl 12 = 0.312±°;°^ [63] follows h' 3 /h 3 « 0.95±g;g||. At the same time this 
Yukawa matrix allows to satisfy severe constraints from lepton-flavour violation, see next 
subsection. We write the Yukawa couplings as hi = \hi\e tlpl including the phases 
Neutrino masses are given in terms of the model parameters as follows: 



|/k» + ^AI~S \h\\\- v/A ^ with A = 2 hifl I ( M i 
K^ + ^Ail- 2 , I/13A3I-— 3— , with A J = ^^-/^- 



M 2 



(6) 



where Am 2 tm = 2.50 x 1CT 3 eV 2 and Am 2 ol = 7.59 x 1CT 5 eV 2 correspond to the squared- 
differences of the eigenvalues of the neutrino mass matrix (3), and the mass difference of 
N\ and N2 is neglected. The third mass eigenvalue is zero due to the flavour structure 
Eq. (5). From Eq. (6) we can estimate the required sizes for the couplings hi and A 5 . 
Assuming I(x) ~ 1 we obtain 



A 5 h 2 Mi VA 



711 



2 



10- 11 (0°) \ 0.05eV 



(7) 



Since hi cannot be too small because of the DM relic abundance, typically A5 has to be tiny 
in order to obtain correct neutrino masses. As discuss later, we impose the perturbativity 
condition \hi\ < 1.5 for the Yukawa couplings. 

2.2 Lepton flavour violation 

Further constraints are imposed on the parameters by limits on charged lepton flavour 
violation. The branching ratios for lepton flavour violating processes i a — > are given 

as 



2 

Br(£ a ^ 7 ) l]n - 



64nG 2 F M* 



t h ^ F2 (w 2 

i=l \ v 



Br (4 -> IpVujJfi) , (8) 



2 If recent indications [62-64] for a non-zero value of the mixing angle #13 should be confirmed 
corrections to Eq. (5) will be necessary. This will be discussed in the last part of Section 2. 
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Figure 1: Contours of Br(^ -» ej) = 2.4 x 10~ 12 in the (M 3 , M n ) plane for various choices of \h 3 \. The 
region to the left of each contour is excluded by /i — > ey. 



where a era = e 2 /(4ir) is the electromagnetic fine structure constant, Gp is the Fermi 
constant and M v is the mass of i] + which we assume to be degenerate with 77° for simplicity. 
The function F^{x) is given by 

, ^ 1 — Qx + 3x 2 + 2x 3 — Qx 2 log x . , 

Fa(x) = 6(T^ • (9) 

The flavour structure of Eq. (5) leads to relaxed constraints from lepton flavour violation 
processes such as /i — > and r — > fij. Because of the two zero's in (5) it follows from 
Eq. (8) that only the third right handed neutrino mass M 3 and the Yukawa coupling h 3 
contribute to // — > e^f process. As a result, r — > /ry gives a more stringent constraint 
than fi — y ey for the neutrino Yukawa couplings h\ , /i2 and the DM mass M\ , and we can 
benefit from the fact that the experimental upper bound Br(r — > fi'-f) < 4.5 x 1CT 8 [66] is 
much looser than Br(/j — > ej) < 2.4 x 1CT 12 [67]. Contours of Br(/i — > ej) = 2.4 x 10 -12 
are shown for several \h 3 \ values in Fig. 1. We take M 3 = 6000 GeV and |/i 3 | = 0.3 as a 
benchmark point in the following discussion. As clear from the figure, for this choice all 
values of M v are allowed, and for M v < 1 TeV we predict \i — > e^y close to the present 
bound. Eq. (7) implies then A5 ~ 10~ 9 . 

Thanks to the restrictions of neutrino oscillation data and lepton-flavour violation 
there are very few independent parameters left. We can choose the following set of four 
independent parameters: 

M„, M x , 6 = M 2 -M 1 , £ = Im(/i^/n) , (10) 

with 8 Mi. Since we fix /i 3 and M 3 to the benchmark point above in order to satisfy 
H — > ej, the relations Eq. (6) determine A 5 as well as \{h\ + h 2 ,)] for a given choice of M v 
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and Mi. However, there is still an undetermined relative phase between hi and h 2l and 
we define the parameter £, which will play an important role in the following. 



2.3 DM relic abundance 

We assume that the lightest right handed neutrino Ni is the lightest of the Z2-odd 
particles, and hence it will be stable and serve as the DM candidate. We assume that 
it is almost degenerated with the second lightest right handed neutrino N 2 . The mass 
degeneracy could be provided by imposing a symmetry for the right handed neutrinos 
Ni. This could be for example the conservation of particle number such that N± and N 2 
form a pseudo-Dirac particle [68]. The smallness of the mass splitting is then related 
to suppressed operators violating the symmetry. The N couple to the SM only via the 
Yukawa interaction with the lepton doublet and therefore our DM is leptophilic. 3 The 
relic density and indirect detection of DM in the model have been investigated with the 
flavour structure of Eq. (5) in refs. [21,22,25,26]. Here we investigate the prospects for 
direct detection of DM in this setup for the first time. 

For the thermal production of DM in this model co-annihilations between Ni and N 2 
have to be considered, since they are assumed to be highly degenerate, leading to an 
enhanced effective annihilation cross section [75]. The effective annihilation cross section 
is written as a e gv = a c s + b e sv 2 + 0(v 4 ). Then the approximate analytic solution of the 
Boltzmann equation which describes the evolution of the DM density is given by 



1.07 x 10 9 x f [GeV- 1 ] . , Mi 

VLh 2 ~ ^- — with x f = — , (11) 

A7*m p i (a e ff + 3b eS / x f ) T f 



where g* is the number of relativistic degrees of freedom at the time of freeze-out Tf 
and m p i = 1.2 x 10 19 GeV. WMAP data [6] implies Qh 2 = 0.11260 ± 0.0036. Taking 
into account co-annihilations of Ni and N 2 we find for the coefficients a e g and b e s in the 
effective annihilation cross section 

e mi , . 

a cS = x - 5 , 12 

^{M 2 + M 2 ) 2 ' 1 ' 

\h 2 + h 2 \ 2 M 2 (M* + M?) e M 2 (M% - 3M 2 M 2 - Mf) 



24tt (M 2 + M 2 Y 2tt (M 2 + M 2 )' 



where the effect of the mass difference between Ni and N 2 is assumed to be negligible. The 
terms proportional to £ 2 come from the co-annihilation process NiN 2 — > £ a ^0, whereas the 
NiNi and N 2 N 2 annihilations lead to the terms proportional to h\ and h 2 , respectively. 



3 Another motivation for leptophilic DM may come from cosmic ray observations from the 
PAMELA [69] and Fermi-LAT [70, 71] experiments, finding an excess of positrons but anti-protons in 
agreement with expectations. In order to obtain the required count rates, however, the annihilation cross 
section must be boosted by a mechanism such as Sommerfeld [72] or Breit-Wigner enhancement [73,74], 
beyond the model considered here. 
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Figure 2: Region in the space of DM mass M\ and £ = |/ii/i2| sin(<^i — ^2) consistent with neutrino 
masses and mixing, lepton flavour violation, perturbativity, and the relic density of DM. The regions 
with different color shadings denoted by A, B, C, D, correspond to different assumptions on M v , with 
A: 2.0 < M7/M1 < 9.8, B: 1.2 < Mr, /My < 2.0, C: 1.05 < M v /M x < 1.20, D: 1.0 < M v /Mi < 1.05. The 
curves show the upper bound on sin #13 from fi — > when the Yukawa matrix Eq. (5) is extended to 
Eq. (14). 

We observe from Eq. (12) and (13) that the s-wave (a e ff-term) is only present due to 
co-annihilations. If there is no phase difference between hi and h 2 , the combination of 
the neutrino Yukawa couplings £ vanishes and only p-wave annihilation remains. This 
corresponds to the helicity suppression for a Majorana fermion. Thus co-annihilations 
and a non-zero phase difference play an important role in obtaining the correct DM relic 
density. 

For the following results we use the micrOMEGAs package [50] to calculate numerically 
the relic abundance of DM. In addition to N\ — N2 co-annihilations, also co-annihilations 
with rj are important, if M v becomes close to Mi. The allowed parameter region in 
the plane of DM mass and the Yukawa coupling £, which is consistent with neutrino 
masses and mixings, lepton flavour violation, and DM relic density is shown in Fig. 2. 
The allowed region is colored and divided into four regions A, B, C, D, corresponding to 
different assumptions on the ratio Mi/M v . The upper bound on £ is imposed by requiring 
perturbativity of the Yukawa couplings. The lower bound on Mi in regions A and B is 
determined by the limit on r — > /ry together with the relic abundance requirement. There 
is no allowed parameter region if M v /Mi > 9.8 because taking into account perturbativity 
as well as r — > /ry the annihilation cross section is suppressed by M*. If M v /M\ comes 
close to 9.8 we are driven to the left-upper corner of the allowed region in Fig. 2. In the 
parameter region C and D we have M^/Mi < 1.2 and co-annihilations with rj become 
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important. Without co-annihilations with rj, the parameter space C and D would not 
appear, and we would obtain a lower bound on |£|. However, if A^i — 77 co-annihilations 
are relevant the correct relic density can be obtained even for vanishing £. In all cases we 
can conclude that the correct relic density is always obtained thanks to co-annihilations 
with either N 2 or 77. 

If M n ~ Mi one may worry about a long-lived charged particle contained in the doublet 
V = {v + >V°)- F° r instance, the predictions for Big Bang Nucleosynthesis (BBN) may be 
altered by the energy injection due to the decay of r/ + into charged leptons [76,77]. We 
have checked that for the parameter ranges of interest 77 decays much faster than 0.01 s 
unless it is degenerate with Ni at the level of 10~ 10 , and hence BBN will be not affected. 

Let us now consider deviations from the flavour structure assumed in Eq. (5). This 
will become necessary if hints for a non-zero #13 [62,64] should be confirmed [65]. In this 
case we expect additional constraints from \x — > ey. Let us consider a small perturbation 
of Eq. (5) as 




h ai = h h 2 h 3 + 0(e 2 ), (14) 



with £j <C hj. Then we allow a non-zero value sin 613 = €3 and deviations of 623 from 7r/4 
as sin #23 = 1/V2 + £4, with e 3j4 <C 1. Diagonalizing the neutrino mass matirx Eq. (3) we 
obtain at linear order in 

1 tan 012^3 nc .s 



fh + e h - x/2(V I h 2 \ ^i + ^)Ai-sec 2 12 ^A 3 



(16) 



If we assume that ei, e 2 and e 3 are real, we obtain from Eq. (8) the following expression 
for \i — > ey. 



Br (jj, — > ey) 



64ttG 2 f M* 



Pt3F2 (^) + V2U n9 MF 2 (^ 



(17) 



Thus a non-zero #13 directly gives a contribution to fi — > ey and €3 = sin 13 is constrained 
by the limit on this process. Using Eq. (6) the parameter P is approximately obtained 
as P w \f2(h\ + hi), and we can obtain an upper bound on e 3 from \i — > ey at each 
point in Fig. 2. 4 Contours of the upper bound on sin #13 are shown in Fig. 2. The upper 
bound becomes severe for small DM mass. Recent results of a non-zero 613 [65] imply 
sin > 0.1 at 3cr. According to Fig. 2 this requires DM masses around the TeV scale 
with f ~ 0(0.1 - 1). 



4 In general, the phase of P depends on the phases of the Yukawa couplings hi and /12, i.e., (pi and cf2, 
not only the phase difference tpi — (p%. For simplicity we set the overall phase of P to zero. This phase 
might play a role if the two terms in Eq. (17) are of comparable size. 



8 



In addition to the extension of the Yukawa matrix Eq. (14) we checked also the effect 
of changing the rl and r2 components into hi +e\ and h 2 + e2- The factor P only changes 
to P ~ 13\/2(/ii + ^l)/6- Moreover, changing the r3 component of Eq. (14) into — (/13 + e), 
the deviation e is required to be zero up to 0(e) from the diagonalization condition of the 
neutrino mass matrix. Therefore these extensions of the Yukawa matrix do not change 
the analysis drastically. 



3 Direct Detection of Leptophilic DM 

3.1 Inelastic Scattering Cross Section 

Inelastic scattering occurs through the effective interactions with quarks which come 
from the 1-loop diagrams shown in Fig. 3. The 3-point vertex effective interactions of N\, 
N 2 and 7 which give a dominant contribution to the inelastic scattering are written as 

C eS = xa^YN^F^ + i (^) W^NxF^ + ic^W^N^ (18) 

where the factor i is a conventional factor to obtain real couplings 012, cu and //12, and 
is the electromagnetic field strength. The coefficient /112 is known as the transition 
magnetic moment between N\ and N 2 . Elastic scattering does not occur through the 
effective interactions because the operators Nij^Nx and Nx<J^ u Nx are identical zero for 
Majorana fermions. General inelastic scattering of DM has been discussed in refs. [68,78], 
and inelastic scattering due to the magnetic moment interactions in ref. [57]. Loop induced 
DM-nucleus scattering for leptophilic DM has been pointed out in ref. [52], and the model 
considered here is a specific realization of "flavoured" DM discussed in ref. [53], where 
similar diagrams to the ones from Fig. 3 have been considered. For another recent model 
for magnetic inelastic DM see ref. [79]. 

In the model considered here, the coefficients ax2, C12 and fix2 are calculated as 

^ Imfeye (Ml m 2 \ 
ai2 ~ Z_, 2 (4tt) 2 M 2 a VM2'M2y<' liyj 

_ v Imfe/ ta i)e (Ml m 2 \ 

^--2^ 2(4tt) 2 M 2 2Ml/m Uf ' Ml I ' (20) 



_ ^ lm(h* a2 h al )e n2r (Ml m|\ 
Cl2 " 2^ 2(4tt) 2 M 2 q c lM 2 'M 2 J' 1 ] 

where g 2 is the momentum transfer and the explicit forms of the function 7 a (x, y), I m (x, y) 
and I c (x,y), which come from the loop integrals, are given in Appendix A. Eq. (5) im- 
plies that Im {h* a2 h a i) = £ and therefore the parameter £ responsible for Nx — N 2 co- 
annihilations controls also the effective interactions of DM with nuclei. 

From the effective interactions, we can obtain three types of differential scattering 
cross sections with a nucleus which has atomic number Z, mass number A, mass rriA, spin 
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Figure 3: Diagrams for the inelastic scattering process of Ni off a nucleus A. 



J a and magnetic moment /i^, see e.g., [53,57]: 



da C c 


Z 2 b\ 2 niA 


dE R 


2nv 2 


da DC 


Zj a cm /i 12 


dE R 


Er 


da DD 




dE R 


7CV 2 



1 - 



E 



R 



1 1 

2m A + Mi 



v 2 Udm v 2 2m A E R 



F 2 {E 



R 



Ja + 1 
3Ja 



F 2 (E 



R 



with the coefficient 



'12 



(ai2 + c 12 q )e. 



(22) 
(23) 
(24) 

(25) 



The cross sections Eq. (22), (23) and (24) are called charge-charge (CC), dipole-charge 
(DC), and dipole-dipole (DD) coupling, respectively. Here E R is the recoil energy, the 
parameter 5 is the mass difference between N 2 and N\ i.e., 5 = M 2 — Mi and /idm is 
the DM-nucleus reduced mass. Magnetic moments of several nuclei are shown in Tab. 1. 
F(E R ) is the nuclear form factor for which we use the parametrization 



F{E_ 



R 



3 [sin(Kr) — nr cos(kx)] 



-kV/2 



(26) 



with « = ^/2m A E R , r = ^R 2 - 5s 2 , R ~ 1.2A 1/3 fm and s ~ 1 fin. F D (E R ) is the 
nuclear magnetic form factor and it is not well-known, see e.g., the discussion in ref. [57]. 
We adopt the following approximation for Fd(E r ). The magnetic moment of a nucleus 
receives contributions from the spin (S n , P ) as well as orbital momentum (L n>p ) of the 
neutrons and protons: 

HA = g s p (S p ) + g s n (S n ) + g\ (L p ) + g{ (L n ) . (27) 
We approximate the magnetic form factor by neglecting the orbital momentum contribu- 
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??Na 




127t 
53 1 


54 ^ve 


133p„ 
55 os 




Ja 


1/2 


3/2 


9/2 


5/2 


3/2 


7/2 


1/2 


Va/hn 


2.629 


2.218 


-0.879 


2.813 


0.692 


2.582 


0.118 



Table 1: Magnetic moments for several nuclei in units of /ijv where hn = e/2m p is the nuclear magne- 
ton [82]. 



tion 5 and use the spin from factors weighted by the corresponding g s factors: 

Fd(Er) ~ g-m^m ■ 1 ' 

The spin-dependent form factors and g s factors are taken from refs. [80,81]. 

In addition to the CC, DC, DD interactions from Eqs. (22), (23), (24) also a charge- 
dipole coupling exists. However there is an additional suppression factor of q 2 compared 
to the other couplings, thus it can be neglected. The DC coupling is singular at Er = 0. 
Therefore the predicted event rate of the DC coupling is enhanced at low recoil energies 
due to the singularity, and we cannot define a total cross section at the zero momentum 
transfer limit ctq C . This situation is the same as in Coulomb scattering. 

3.2 Comparison of the Predicted Event Rate with Experiments 

We compare the event rate calculated from the effective interactions with XENON100 
[42], KIMS [48] and DAMA [44] data. The DD coupling might be important for KIMS or 
DAMA [57] since in these experiments, the target nuclei are iodine (I) and cesium (Cs) for 
KIMS, iodine and sodium (Na) for DAMA, which have a large nuclear magnetic moment 
as can be seen from Tab. 1. The event rate is written as 

dR v-^ 1 f da „. . o . . 

wr^M^L. M-/ f{v)d v ' (29) 

U nuclei 1 Qct Jv>v ™™ « 

where p & ~ 0.3 GeVcm -3 is the local DM density, Md e t is the mass of target material, 
v min is the minimum velocity required for DM to scatter off a nucleus with recoil energy 

^ = -^=(^ + 5), (30) 



y / 2m A E R V Pdm 

and f(v) is the local DM velocity distribution function in the rest frame of the Earth. 
It is obtained by a Galilean transformation from a Maxwell-Boltzmann distribution in 



5 The ratio of spin and orbital contributions to the magnetic moment in Eq. (27) are 0.59 : 0.41 for 
Sodium, 0.52 : 0.48 for Iodine, 0.96 : 0.04 for Xenon, -0.38 : 1.38 for Cesium. Therefore, neglecting 
the orbital contribution is an excellent approximation for Xenon. For the other nuclei this introduces an 
error of about an factor 2 and therefore the limits derived from KIMS and DAMA should be considered 
only approximate. 
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Energy range 


Quenching factor 


XENON100 


8.4 - 44.6 keV 




KIMS 


3.6 - 5.8 keVee 


0.1 (Cs), 0.1 (I) 


DAMA 


2-8 keVee 


0.3 (Na), 0.09 (I) 



Table 2: The energy range and the quenching factor for the experiments XENON100 [42], KIMS [48], 
and DAMA [44]. For XENON100 we use the same light-yield function L cS as in ref. [42]. 



the rest frame of the galaxy with the velocity dispersion vq = 220 km/s and the escape 
velocity from the galaxy v csc = 544 km/s. The velocity distribution function f(v) is 
normalized to J f(v)d 3 v = 1. The relative velocity of the Earth to the galaxy is v e = 
% + f rbCOS7cos [27v(t — i )/year] with v = v + 12 km/s, t> or b = 30 km/s, cos 7 = 0.51 
and to = June 2nd. We must evaluate the following velocity integrals to predict the event 
rate: 

[<*> f{v+Ve) , 3 

ClKnin,We) = / <TV, (31) 

Jv ■ V 
/>oo 

(2(v min ,v e ) = vf(v + v e )d 3 v. (32) 

" ^min 

The analytic formulas for the DM velocity integrals given in refs. [60, 83] are used. The 
total predicted event rate in the XENON100, DAMA, and KIMS experiments is obtained 
by integrating the differential event rate with respect to an appropriate recoil energy range. 
We use the energy range and the quenching factors shown in Tab. 2. The quenching factor 
is the ratio of the energy deposited in scintillation light to the total nuclear recoil energy. 

In Fig. 4 we illustrate the relative importance of the the CC, DC, DD interactions from 
Eqs. (22), (23), (24) for the XENON100, KIMS, and DAMA experiments by calculating 
the total event rate induced from each of the three interaction types separately. We 
observe from the left panel that typically CC interactions are more important for small 
masses Mi, which follows from the different dependence on the DM mass of 612 and [L\i- 
The value Mi where CC becomes subdominant depends on the ratio M n /M\. The right 
panel of Fig. 4 shows that for XENON100 the DC coupling is more important, whereas 
for KIMS and DAMA DD dominates, because of the large magnetic moments of iodine 
and sodium. The features of the DAMA curves around Mi ~ 20 GeV in both panels are a 
consequence of the presence of the two elements (I and Na) with rather different masses. 
In general the relative importance of CC, DC, DD depends on the ratio M v /M\ and to a 
lesser extent on 5. The main conclusion is that depending on the region in the parameter 
space and depending on the considered experiment, any of the three interaction types can 
be important and all of them have to be taken into account. 

In order to derive constraints on the model we calculate the total event rate for 
XENON100 and KIMS in the energy range given in Tab. 2 and require that the pre- 
dicted rate is less than 0.0017, 0.0098 kg'May" 1 for XENON100 and KIMS, respectively. 
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Figure 4: Relative contributions of the charge-charge (CC), dipole-charge (DC), and dipole-dipole (DD) 
interactions to the total predicted event rate in XENON100, KIMS, and DAMA. The left panel shows 
the contribution from CC relative to the sum of DC and DD, the right panel shows the ratio of the DC 
and DD contributions. We assume M v /Mi = 1.5 and S = 0. 



The upper bounds are obtained from the observed 3 events with 3a of the statistical 
error in the 48 kg fiducial volume during 100.9 live days exposure in the signal region 
for XENON100 [42], and from ref. [48] for KIMS. For DAMA we perform a x 2 fit to the 
modulation amplitude in bins of observed scintillation energy between 2 and 8 keVee. 
In Fig. 5 and Fig. 6 we show the bounds from XENON100, KIMS and allowed regions 
from DAMA for the coefficients &i 2 (see Eq. (22)) and fi± 2 (see Eqs. (23), (24)), respec- 
tively. These bounds are compared to the regions as predicted in the model according to 
Eqs. (19), (21), (25) for b\ 2 and Eq. (20) for /ii 2 . The colored regions correspond to the 
regions shown in Fig. 2, satisfying constraints from neutrino masses and mixing, charged 
lepton-flavour violation, the relic DM density, and perturbativity. The ratio M v /Mi is 
taken in the range 1 < M v /Mi < 9.8, with the same color shading as in Fig. 2. There is 
no allowed parameter space for M^/Mi > 9.8, as discussed earlier. 

We observe that the values of | &12 1 and \fJ*i 2 \ obtained in this model are too small 
to account for the signal in DAMA. For very small mass splittings 5 between Aq and 
N 2 some regions of the parameter space are excluded by XENON100 data. The con- 
straints become weaker for larger S, since increasing inelasticity suppresses the scattering 
event rate. Relatively large values of | &12 1 are obtained for close to degenerate i\q and 
77, Mrj/Mi < 1.05 (dark-red region), because of the behavior of the function I a (x,y) near 
x = 1, where I a (x, y) ~ y^ 1 and y = m 2 a /M^ is small. The region excluded by XENON100 
for Mi ~ M v ~ 2 TeV becomes allowed for 5 > 120 keV (bottom- right panel). By com- 
paring Figs. 5 and 6 we observe that the model predicts values of [/Xxa | too small to be 
tested by current direct detection data. The enhancement for the transition magnetic 
moment I//12I for M n /Mi < 1.05 is less than for |6 12 | due to a different behavior of the 
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Figure 5: Bounds from XENON100, KIMS and allowed regions for DAMA in the (Mi,|6i 2 |) plane 
(charge-charge interaction). The mass difference 6 is taken as keV (the left top panel), 40 keV (the 
right top panel), 80 keV (the left bottom panel) and 120 keV (the right bottom panel). The shaded 
regions correspond to the values of 612 predicted in the allowed parameter space of the model, as shown 
in Fig. 2, with the same color shading for different values of the ratio M v /Mi. 



loop functions I a (x,y) and I m (x,y). 

We conclude that current data from XENON100 start to exclude some parameter space 
of the model, in case of degenerate configurations Mi ~ M 2 ~ M v ~ few TeV. In Fig. 7 we 
show the regions excluded from XENON100 overlayed to the globally allowed regions from 
Fig. 2 as dark blue, by translating the the | &12 1 constraint into a bound on £. Furthermore 
we show in Fig. 7 the estimated sensitivity for XENON1T. Using the sensitivity for the 
elastic WIMP-nucleon scattering cross section from ref. [84] we estimate that XENON1T 
will constrain the event rate to be less than 1.59 x 10~ 5 [kg _1 day -1 ]. (We assume the 
same nuclear recoil energy range as for XENON100.) Then we compare this number to 
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Figure 6: Bounds from XENON100, KIMS and allowed regions for DAMA in the (Mi, \fi 12 \) plane 
(dipole-charge and dipole-dipole interaction). The mass difference <5 is taken as keV (the left top 
panel), 40 keV (the right top panel), 80 keV (the left bottom panel) and 120 keV (the right bottom 
panel). The shaded regions correspond to the values of /X12 predicted in the allowed parameter space of 
the model, as shown in Fig. 2, with the same color shading for different values of the ratio M v /Mi. 



the event rate induced in the model assuming several values for the mass splitting 5. From 
Fig. 7 we find that for 5 < 40 keV future data from the XENON1T experiment [84] will 
dig deeply into the allowed parameter region of the model. For 40 keV < 5 < 120 keV the 
degenerate region Mi ~ M 2 ~ M v ~ few TeV will be tested. We note however, that no 
signal is guaranteed for direct detection. In the N% — r\ co-annihilation region (dark- and 
light-red regions, where M n jM\ < 1.2) no lower bound on the parameter |£| is obtained, 
leading to arbitrarily small values of \bi2\ and \fin\, which implies a vanishing signal in 
direct detection experiments. 
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Figure 7: Same as Fig. 2 with constraints from XENON100 (blue) and sensitivity from XENON1T 
(gray) overlayed. We assume 6 =0 (left top), 40 keV (right top), 80 keV (left bottom), 120 keV (right 
bottom). 



3.3 Monochromatic Photon from the Decay of N2 

The excited DM state N 2 decays to Ni and a photon through the transition magnetic 
moment. The diagrams of the decay process are shown in Fig. 8 and the decay width is 
calculated as 

T(N 2 N a ) = ^5 3 . (33) 

7T 

Notice that the effective interaction b\2 does not contribute to the decay width since the 
emitted photon is on-shell. The decay of N 2 produces a monochromatic photon of energy 
E 1 ~ 5. If the decay happens inside a DM detector this monochromatic photon would 
contribute to the electromagnetic event rate. Although typically such events are rejected 
in order to search for nuclear recoils one may be able to place constraints on the model 
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by requiring that the electromagnetic event rate induced by the decay of N 2 has to be 
less than the observed rate. A similar mechanism has been used in ref. [85] in order to 
explain the DAMA modulation signal. 

Following ref. [85] we estimate the photon induced event rate in the model under 
consideration for the XENON100 experiment. The excited state N 2 is produced by the 
inelastic scattering with nuclei inside the Earth which is composed of various elements 
such as Fe, O and Si. The event rate in XENON100 is given by 

§^ = 17^- £ / d*v^vf(v) [ <? rni (r)P(r,v), (34) 

dh R Mipxe i= ^ lei Jv>v min d^R ./Earth 

where pxe is the mass density of the XENON detector 21.9 g/cm 3 , <j\ is the total inelastic 
scattering cross section which includes the charge-charge, dipole-charge and dipole-dipole 
interactions, and rii(r) is the number density for the given nucleus i inside the Earth. 
Note that E R is the nuclear recoil energy in the N\ + A — > N 2 + A scattering process. 
The contribution of the dipole-dipole interaction is much smaller than the ones from the 
charge-charge and dipole-charge interactions since the fraction of isotopes with a sizable 
magnetic moment in the Earth is less than a few %. In Eq. (34), P(r, v) is the probability 
that an N 2 which is produced by the scattering of DM with velocity v at the position r 
decays inside the XENON100 detector. It is given by 

P(r,v) = — -! e -r|r-r*l/«, (35) 

4n(r - r Xc ) z Vf 

where vj = ^/v 2 — 2(5 + Ep)/Mi is the velocity of the produced N 2 , and rx e is the 
position of XENON detector on the Earth. The total gamma event rate i? 7 is obtained 
by integrating Eq. (34) over the recoil energy E R . 

In order to obtain a rough estimate of the induced event rate we introduce some 
approximations. We use the averaged number density of the elements in the Earth n ~ 
9.85 x 10 22 cm" 3 , the averaged atomic number Z ~ 29.9 and the averaged magnetic 
moment of nuclei JL A / '(jln — 3.46 x 1CT 2 , which are calculated by taking into account the 
structure of the Earth such as the crust, mantle and core [85]. Replacing n^(r) by its 
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average, it can be pulled out of the r-integral in Eq. (34) and the integration is performed 
analytically: 



'Earth 
16 



d 3 rP(r,v) = - 



( e -*rr<B/v, _ 1) + 1 
21> ffi K ' 



(36) 



where r e = 6.4 x 10 m is the radius of the Earth. The remaining integrations over v and 
Er are done numerically. 

With this approximation we estimate the total predicted event rate in XENON100 
for typical parameters of the model. We find that the maximal rate is approximately 
R m ax _ 2.0 x 10~ 7 kg _1 day _ \ when 8 ~ 40 keV and M^jM^ = 1. This result should 
be compared with 22 events obtained in the electromagnetic band in the 40 kg fiducial 
volume during 11.17 live days exposure in the DM search window by XENON100 [41]. 
Hence, since the predicted event rate is several orders of magnitude smaller we conclude 
that the monochromatic photon from the N 2 decay will not lead to any observable signal 
in DM direct detection experiments. 



4 Summary and Conclusions 

We have considered a model proposed by Ma [10], providing an economical extension 
of the Standard Model to accommodate neutrino masses and DM. The Standard Model 
is extended by three fermion singlets Ni ("right handed neutrinos") and an inert scalar 
doublet i], where the new particles transform odd under a Z2 symmetry, making the 
lightest of them a stable DM candidate. In our case N% is the DM particle. We investigate 
the parameter space of the model consistent with neutrino masses and mixings, bounds on 
charged lepton-flavour violation, perturbativity, and the correct relic DM abundance due 
to the thermal freeze-out mechanism. We find that in order to obtain the correct relic DM 
abundance co-annihilations are always important, either between the two lightest fermion 
singlets Ni and N 2 , or between Ni and the inert doublet r\. 

In this model DM has no direct couplings to quarks and gluons. Despite this leptophilic 
nature of DM, scattering off nuclei is possible at 1-loop level by photon exchange. We 
have calculated the relevant loop processes in an effective field theory approach. One 
obtains effective charge- charge, dipole-charge, and dipole-dipole interactions between DM 
and nuclei, leading to a non-vanishing scattering rate in DM direct detection experiments. 
The scattering is inelastic and in order to obtain a sizable scattering rate N\ and N 2 have 
to be highly degenerate, with mass differences 8 less than few 100 keV. This is consistent 
with the need for co-annihilations to obtain the correct relic abundance. Although the 
scattering cross section in this model is too small to account for the DAMA annual 
modulation signal, we find that for mass differences 8 < 120 keV current data from 
the XENON100 experiment start to exclude certain regions of the parameter space. The 
predicted event rate for XENON100 is dominated by the charge-charge interaction. Future 
data, for example from XENON1T, will significantly dig into the allowed parameter space 
and provide a stringent test for the model provided 8 is small enough. 
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Note added. 

After this work has been completed the Daya Bay reactor experiment released their 
data [65], establishing a non-zero value of #13 at more than 5a with sin 2 2^13 = 0.092 ± 
0.016 ± 0.005. For non-zero values of #13 additional contributions to fi — > cy are induced, 
providing further constraints on the model. Daya Bay data imply sin# 13 > 0.1 at 3a, 
which constrains DM masses around the TeV scale, see Fig. 2. 
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Appendix A 

Explicit Functions for the Effective Interactions 

Here we give the explicit functions for the effective interactions. The functions I a (x, y) 
and I m (x) are given as follows, 



3u 2 - 6u + 1 



3 J xv? — (1 + x — y)u + 1 ' 

I m ( X ,y) = -f . 

J xu 2 — (1 + x — y)u + 1 

The analytic formulas of these integrations are given as follows. 

(i) If (1 + x - yf - Ax > 0, 



(37) 
(38) 



1 

x 
1 

X 



1 + 



4 (x, y) 

lm(x,y) 

(ii) If (l + x-?/) 2 -4x = 

1 



3Al - QA+ + 1 



3{A + - A_) 



log 



A+-1 



A, 



3A 2 _ - 6A_ + 1 



3{A + -A_) 



log 



A_ - 1 



A_ 



A,-l 



A. 



A, - A_ 



log 



A_ - 1 



A_ 



,(39) 
(40) 



h (x, y) 
Im(x,y) 
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1 

X 



1 + 2(A - 1) log 

2 + (2A - 1) log 



A -l 



A Q 
A. -I 



+ 



3A 2 - QA + 1 
3A (A - 1) 



A, 



(41) 
(42) 
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(iii) If (1 + x - yf - 4a; < 0, 



40, y) 



Im{x, y) 



1 

x 



1 

X 



\ , B + + B_-2 ^ 
1 H o log 



6(5+ - - 1) -4 



Bl + B 2 _ 



3(5+ - B_ 
1 H log 



Tan 



-i 



B + - B_ 



[B A 



B\ + B 2 - B+- B_ 
l) 2 + (5_-l) 2 



(43) 



B\ + B 2 _ 



(2B+ - 1)(2B_ - 1) - 1 



2(B + -B_) 



Tan 



-i 



B+-B. 



Bi + B 2 _-B+- B_ 



• (44) 



A±, A Q and B± are denned as 

A ± 
A 
B ± 



1 + x — y ± a/(1 + s — y) 2 — 4x 



2x 

1 +x — y 

2x ' 

1 + x — y ± a/4x — (1 + x — 7/)^ 
2x 



The function I c (x,y) is the same as I m (x,y). These functions are continuous and smooth 
for < x, y < 1. For ~ y <C x <C 1, these functions approach to 



1 2 

iafoy) ->■ - + - logy, 



Im(x,y) ->• 



1 



(45) 
(46) 



Therefore, the obtained parameters | &12 1 and at lowest order agree with the result of 
ref. [53] where the parameter A 2 in ref. [53] corresponds to Im(/i* 2 /i Q i)/2 in our notation. 
The difference of the relative sign comes from the denifition of the effective operators. 
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